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Abstract 

Holographic functional methods are introduced as probes of discrete time-stepped maps that lead 
to chaotic behavior. The methods provide continuous time interpolation between the time steps, 
thereby revealing the maps to be quasi-Hamiltonian systems underlain by novel potentials that 
govern the motion of a perceived point particle. Between turning points, the particle is strictly 
driven by Hamiltonian dynamics, but at each encounter with a turning point the potential changes 
abruptly, loosely analogous to the switchbacks on a mountain road. A sequence of successively 
deepening switchback potentials explains, in physical terms, the frequency cascade and trajectory 
folding that occur on the particular route to chaos revealed by the logistic map. 



1. INTRODUCTION 



In a previous paper [1| we have discussed how functions of position, defined on a discrete 
lattice of time points, may be analytically interpolated in time by functions of both x and t, 
defined on a continuum of time points, framed by the lattice, through the use of solutions to 
Schroder's nonlinear functional equation [2]. For a more pedagogical discussion to augment 
this Introduction, we encourage readers to look at that companion paper. 

If the effect of the first discrete time step is given as the map x i— y fi (x,s), for some 
parameter s, Schroder's functional equation is 

s*(x,s) = *(f 1 (x,s),s) , (1) 

with \I/ to be determined. So, 

f 1 (x,s) = ^- 1 (s^(x,s),s) , (2) 
where the inverse function obeys Poincare's equation, 

y- 1 (sx,s) = f 1 (y- 1 (x,s),s) . (3) 

The nth iterate of ([!]) gives 

s n V(x,s) = *(/i(...A (/x (x,s), s)... ,s),s) , (4) 

with f\ acting n times, and thus the nth order functional composition — the so-called 
"splinter" of the functional equation — 

f n (x, s) = A (- • • h (x, s) , s) ■ ■ ■ , s) = M/- 1 ( x , s ) jS ) . (5) 

A continuous interpolation between the integer lattice of time points is then, for any t, 

f t (x,s) = ^- 1 (s t ^(x,s),s) . (6) 

This can be a well-behaved and single-valued function of x and t provided that ^l/ -1 (x, s) is 
a well-behaved, single- valued function of x, even though \1> (x, s) might be, and typically is, 
multi-valued. 

As discussed in the interpolation can be envisioned as the trajectory of a particle, 

x(t) = f t (x,s) , (7) 

where the particle is moving under the influence of a potential according to Hamiltonian 
dynamics. The velocity of the particle is then found by differentiating ([6]) with respect to t, 

^- = (\ns)s t ^(x)(^- 1 (s t ^(x)))' , (8) 

where any dependence of ^ on s is implicitly understood. Therefore, the velocity will inherit 
and exhibit any multi-valuedness possessed by \l/ (x). 

Indeed, suppose that Hamiltonian dynamics have been specified, trajectories have been 
computed, and fi = x(t = 1) has emerged in terms of initial velocity and initial position 
x (t — 0) = x. A solution of the functional equation ([T]) can then be constructed, and 
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expressed in very physical terms, as just an exponential of the time elapsed along any such 
particle trajectory that passes through position x. That solution is 

*(x) = s T ^ , T(x) = J X -^- ) . (9) 

Here v is the velocity as a function of the position along the trajectory. (Hence v (x (t)), 
if you wish. Of course, for a Hamiltonian system with time translational invariance, we 
may supplant x (t) — >■ x (t + At) in all these expressions, as well as those to follow.) For a 
particle moving in a potential V (x) at fixed energy, E, the velocity can be expressed in the 
usual way as v (x) = ±y/~E — V (x), with suitably chosen mass units. 

When written in a form more closely related to Schroder's functional equation, the solu- 
tion (jUJ) is simply 

m (f t ( x )) = S J * * (x) . (10) 



At t = 1, f£ '^yj = I; an d Schroder's equation (JTj) re-emerges. But, as this construction 
clearly shows, one must be careful at turning points where v — 0, especially if these are 
encountered at finite times along the trajectory. Typically these turning points produce 
branch points in \l/ so that it is multi-valued. 

An interesting physical effect of such branch points, when they are encountered in finite 
times, is the possibility to switch from one branch of the underlying analytic potential 
function to another, thereby changing the functional form of the Hamiltonian governing the 
motion along seemingly identical intervals of a real-valued trajectory. We will therefore call 
these gwasz-Hamiltonian systems, and we stress that when such switches occur this is not 
standard textbook Hamiltonian dynamics. Further explanation will be given below in the 
context of the logistic map. 

The interpolation ([6]) can also be viewed as a holographic specification on the x, t plane 
[jjj], determining f t (x) in the surrounded "bulk" from the total data given at the bounding 
times, {x} U {fi (x)}. In this point of view, fixed points in x complete the boundary data 
and facilitate the solution of Schroder's equation through power series in x, hence leading 
to V (x) which was not known a priori [l|. 

In this holographic approach, the potential first appears as a quadratic in ty/ty', 

^) = -(M 2 (|f|)\ (ID 

up to an additive constant, where any dependence of ^ on s is again implicit. The x 
dependence of the potential, 

V{x) = -v 2 {x) , (12) 
follows from that of the velocity profile of the interpolation, defined and given by 

/ x df t (x) 

V (X) EE 



= (M tfCsH*- 1 (*(*)))' = ^Wm- ( 13 ) 



dt 

Monotonic motions between two fixed points, such as occur for the Ricker model [5[, provide 
the most elementary examples [l[. 

However, situations where one or both of the fixed points are absent were not fully 
addressed in our previous work. This is precisely the situation that occurs when turning 
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points are encountered at finite times in the particle dynamics interpretation, and leads to 
an intriguing modification in the physical picture involving the potential and its effect on 
the particle trajectories. We consider here a specific example of such a situation involving 
the well-studied logistic map of chaotic dynamics, namely 



ipie oi ; 



x i — y sx (1 — x) . (14) 

Applying our functional methods to this example, the resulting interpolations from a discrete 
lattice of time points to a time continuum then allow us to appreciate analytic features of 
the logistic map, and to derive the governing differential evolution laws — indeed, subtly 
time-translation-invariant Hamiltonian dynamical laws — of the underlying physical system, 
hence to obtain potentials that were not previously known. 

For the all-familiar logistic map illustrating transition to chaos, the functional interpo- 
lation reveals that there are well-defined expressions for the continuous time evolution of 
this map for all parametric values of the map, whether chaotic or not. We obtain agree- 
ment with the explicit closed-form solutions of ([1]) for the special values of the parameters: 
s = —2, 2, and 4. (These explicit solutions have been known for almost a century and a half 
0.) Moreover, as indicated, from ffTTJ we can now find the potentials needed to produce 
these explicit, continuously evolving trajectories in the language of Hamiltonian dynamics. 

A new feature for the potentials so obtained for the s = 4 case in particular, and also 
for other values of s > 2, is that they must change at discrete intervals of the envisioned 
Hamiltonian particle's motion, to be consistent with the evolution trajectories: Every time 
the particle hits a turning point, the potential changes [TlJ. We therefore call the corre- 
sponding V^s "switchback potentials" and to echo our previous general remarks, we will call 
the dynamical system as a whole "quasi-Hamiltonian." For s = 4 the switchback potentials 
deepen successively and thus lead to more rapid, higher frequency motion, resulting in the 
familiar chaotic behavior of the discrete s = 4 logistic map that they interpolate. The famil- 
iar frequency-cascade-and-folding behavior of the chaotic discrete map is thus understood 
from — indeed, explained by — the subtle properties of the switchback potentials. 



2. THE LOGISTIC MAP 

Consider in detail the logistic map ffT4"j) on the unit interval, x G [0,1]. Schroder's 
equation for this map is 

(x, s) = * (sx (1 - x) , s) . (15) 
The inverse function satisfies the corresponding Poincare equation, 

(sx,s) = s^- 1 (x,s) (l-^- 1 (x,s)) . (16) 

As originally obtained by Schroder, there are three closed-form solutions known, for s = —2, 
2, and 4: 

V?> f r ( 1 \ \ 1 ( 2x ir 



^ (x, -2) = — \2tc - 3 arccos [x , (x, -2) = cos — = + - 

1 ' ; 6 V V 2JJ 2 Vv/3 3 

* (x, 2) = - 1 - In (1 - 2s) , (x, 2) = ~ (l - e~ 2x ) , 

\1/ (x, 4) = (arcsin 

V^) , * _1 (x,4) = (sinv^) • 



4 



Note that, while the \I/ are multi-valued, the inverse functions are all single-valued. 
More generally, consider a power series for any s. 



ijr-i^g) = x + x^2x n c n (s) . (18) 

n=l 

The Poincare equation then leads to a recursion relation for the s-dependent coefficients. 

1 n 

Cn+1 = i 7n+^ y~! c i Cn -i ' (^) 

j=0 

with Co = 1, ci = C2 = tjz^wz^) ? e tc- The explicit coefficients are easily recognized 
for s = —2, 2, and 4, and immediately yield the three closed-form cases. Similarly let 

oo 

V(x,s) = x + x^(-x) n d n {s) . (20) 

n=l 

Then, as a consequence of Schroder's equation, d\ — 1/ (1 — s), and for n > 2, 

where |_' ' " J is the (integer-valued) floor function. In principle, these series solve ( !T6|) and 
( f!5|) for any s, within their radii of convergence. 

From extensive numerical studies we believe the radius of convergence for ( 12"U|) depends 
on s as follows. 

( 1 

1-- if0<s<2, 

= < S (22) 

- if 2 < s < 4 . 
I 4 

Preliminary work further suggests that these conjectural results can actually be established 
analytically by comparison with known convergent series, at least in some cases. 

The radius of convergence for the inverse function series ( TTBl) is intriguing for < s < 1 . A 
fit by Mathematica to radii obtained numerically from explicit series coefficients c n (s) up to 
n = 200, for 20 different values of s in this range, using the functional form | (1 — s) a exp bs, 
gives a = 0.78 and b = 0.62. Thus R^-i (s) falls from 1/4 to zero as s goes from to 1. 
However, for the parameters of most interest, namely 2 < s < 4, numerical studies suggest 
the series all have infinite radius of convergence, indicating the are entire functions for 
these parameters. This is true for s = 2 and s = 4, of course, as is evident from (fT7|) . The 
numerics provide compelling evidence that this is also true for all intermediate s. 

Given convergent series, the functional equations can then be used to continue the series 
solutions outside their radii of convergence and thereby obtain accurate determinations of the 
various branches of and of the large-but-finite argument behavior of \J r_1 , in a manner 
familiar from, say, the T and ( functions. (See [16| for details, and the Appendix for a 
numerical example.) Comparison of the results from these series-plus-functional-equation 
methods with the known exact solutions flT7|) shows perfect agreement. 
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The first few terms for \l/ and \I/ 1 for generic s are given explicitly by [l2| 

. » x 2 2s x 3 s(l + 5s 2 ) x 4 

* (ar, s) = x + + — "~ir + 



s-1 (s + l)( s -l) 2 ^ 2 + 8 + l)( 8 + l)(a-iy 

2s 3 (3 + 2s + 7s 3 ) x 5 

+ (s 2 + l)(s 2 + s + l) (s + l) 2 (s-l) 4 

2s 3 (1 + 3s + 14s 3 + 14s 4 + 7s 5 + 21s 7 ) x 6 ,~ 

+ (s 4 + s 3 + s 2 + s + 1) (s 2 + 1) (s 2 + s + 1) (s + l) 2 (s - l) 5 + ^ ' <2o) 

_ i , . x 2 2 x 3 5 + s x 4 

(x, s) = x H h 



1-s ' (s + l)( s -l) 2 (s + l)(s 2 + s + l)(l- s ) 3 

2 (7 + 3s + 2s 2 ) x 5 
+ (s 2 + l)(s 2 + s + l) ( s + l) 2 ( s -l) 4 

2 (21 + 14s + 14s 2 + 8s 3 + 3s 4 ) ^ + (x 7 ) (24) 

(s 4 + s 3 + s 2 + s + l)(s 2 + l)(s 2 + s + l)(s + l) 2 (l-s) 5 1 ' ' 

from which we infer that $ (x, s) fx and (x, s) /x are actually series in xj (1 — s) with 
s-dependent coefficients that are analytic near s — 1. 

The trajectories interpolating the splinter (integer t) of the logistic map are then, 

^ T w tr / * \ t s'(l-s') 2 2s*(l-s*)(s-s*) o 

X (t) = (s'tf X, S , s) = S* X + — i - X 2 + i 5-^ X 3 

s'(l-s*)(s-s')(l + 5s 2 -(s + 5)s*) 4 



X 



(s + l)(s 2 + s + l) (s-1) 3 
2s* (1 - s*) (s - s*) (s 2 - s r ) (7s 3 + 2s + 3 - s* (2s 2 + 3s + 7)) 5 
(s + l) 2 (s 2 + l)(s 2 + s + l) (s-1) 4 X 



+ 0(x 6 ) . (25) 



The trajectories are single-valued functions of the time so long as \1/ 1 is single-valued, and 
in fact, they exist even for s — > 1 as formal series solutions. Explicitly, 

lim (s'tf (x, s) , s) = x - i x 2 + t (t - 1) x 3 - -t (t - 1) (2t - 3) x 4 

+ ~£ (t - 1) (t - 2) (3t - 4) x 5 - — t (t - 1) (t - 2) (I2t 2 - Alt + 31) x 6 
3 X2 

+— t (t - 1) (t - 2) (30t 3 - lilt 2 + 302i - 157) x 7 + O (x 8 ) . (26) 
30 

For the three special cases, s = — 2, 2, and 4, there are closed-form results for various 
quantities of interest. For example, for s = 4, the trajectory and velocity are given by 

x(t)| s=4 = (4* tf(x,4) ,4) = (sin (2* arcsin v/x")) 2 , (27) 
(as originally presented in fy, p 306) and by 

—x (t) = (2 1+t In 2) sin (2* arcsin ^fx) cos (2* arcsin a/x) arcsin sfx (28a) 

(JjL 

= (In 4) y/x(t) (1 -x(t)) arcsin y/x(t) . (28b) 
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The last expression evinces a continuous time-translational invariance. However, this ve- 
locity function has branch points (i.e. turning points) at x (t) =0 and x (t) = 1, so some 
care is needed to determine which branch of the function is involved, particularly when the 
turning points are encountered at finite times. The system is therefore a quasi-Hamiltonian 
one, as we have previously indicated, and as we shall discuss in more detail. 

With this caveat in mind, we may thus deduce the velocity profile v (x) = 

effective potential, and the force for the underlying quasi-Hamiltonian system: 

dx (t) 



, the 

t=o 



r(.r) = — — — = (In 4) \J x (1 — x) arcsin \fx , (29) 

V (x) = — v 2 (x) = (ln4) 2 x (x — 1) arcsin 2 \fx , (30) 
F (x) = ( x ) = (hi4) 2 (arcsin \/x) (\/ x (1 — x) — {2x — 1) arcsin \/x^j . (31) 



Note x = is an unstable fixed point for the system. Also, the system is time-translationally 
invariant, so these same expressions hold at all times everywhere along a trajectory, provided 
due care is taken to determine which branch of the various functions is in effect. Thus 
V (x (t)) = (ln4) 2 x (i) (x (t) - 1) arcsin 2 ^x (t), etc. 

Similar closed-form results hold for the other two solvable cases, s = —2 and s = 2. In 
principle, there are also potentials underlying the trajectories for other s by dint of the series 
constructions. (We illustrate one other case, s = 3, in the Appendix.) 

From the general expression for x (t) in terms of \l/ and \fr , and with use of the chain 
rule, the potential may be expressed entirely in terms of (ln^)', as given in the Introduction 
by ffTTj) . The Schroder auxiliary function is then recognized as just an exponential of the 
time function 

rw ^7*r (32) 

computed along a zero-energy trajectory, as given by Q). 

3. NOVEL POTENTIALS AND SWITCHBACK EFFECTS 

The effective potentials for all three closed-form cases are somewhat unusual: 

V (x, s = 4) = (In 4) 2 x (x — 1) arcsin 2 y/x , (33) 

V ( x , s = 2) = - (in V2^j 2 (1 - 2x) 2 In 2 (1 - 2x) , (34) 

V(x,s = -2) = (ln(-2)) 2 (2x + 1) (2z - 3) ^27r-3arccos (x - . (35) 

Another way to express the potential for s = 4 is similar in form to that for s = —2, namely, 

V (x, s — 4) = (In 2) x (x — 1) (n — arccos (2x-l)) 2 . (36) 

Indeed, it is well-known that the logistic maps for s = 4 and s = —2 are intimately related 
through the functional conjugacy of the underlying Schroder equations. But note the 
potential for s = —2 is in fact complex, since In (—2) = In 2 + as are the trajectories 
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under real time evolution for this case. Of course, since the complexity of V (x, s — —2) is 
solely a multiplicative factor, if we switch to complex time, r = (In 2 + m) x t, then 



d 2 x = d V(x,s = -2) 

dr 2 8x 2(ln2 + m) 2 1 ' 

again yields real trajectories x(t). As trivial as this observation is, nevertheless it does 
raise several questions about complex x, and about the behavior of V in the complex plane. 
We discuss only one aspect of this behavior here, and consider non-principal values for the 
multi- valued functions v (x,s = 4) and V (x, s = 4). In particular, all branches of the arcsin 
function are important to understand the behavior of the explicit trajectories (|2"7|) . 

Switchbacks on the road to chaos. An interesting new feature, which we shall call 
the switchback effect, appears for a particle moving in the s = 4 effective potential. This 
is a distinguishing feature that we encounter for the s = 4 chaotic map, but not, say, for 
the non-chaotic s = 2 map. The effect is traceable to the fact that, while the trajectory 
is single-valued as a function of the time, the velocity and hence the effective potential are 
not single-valued as functions of position. Moreover, the branch points of the multi-valued 
functions are encountered by zero-energy trajectories at finite times for the s = 4 case, 
unlike, say, for the s = 2 case. For the latter case, it takes an infinite time for a particle 
with zero energy to reach a turning point. 

Switchbacks are essentially transitions from one branch of the position-dependent velocity 
function to another, and occur at the turning points encountered at finite times by the 
interpolating particle trajectories, (|27p . The effect is easily seen upon viewing animations 
of these trajectories fl3j | . 

Consider a particle with zero total energy that starts in the V (x, s — 4) potential given 
above, in the region < x < 1. If the particle is initially moving to the left, it will continue 
towards x = 0, taking an infinite amount of time to reach that turning point. But if the 
particle is initially moving to the right, it will reach the x = 1 turning point in a finite 
amount of time that depends on its initial x, as given by 



. , s 1 f 1 dy 1 / tt/2 \ 

At (x) = - — / 7 = - — In '- — = 

ln 4 Jx y/y(l-y) (arcsin y/y) In 2 \ arcsin y/x J 



(38) 



So, for example, starting from the midpoint x — 1/2, it takes unit time to reach the turning 
point at x = 1. 

Upon reaching x = 1, the explicit form of the time-dependent solution (1271) exhibits the 
classical counterpart of a sudden transition that keeps E = 0, but changes the potential for 
the return trip towards x — 0, exactly as follows from the particle moving on a different 
branch of the arcsin function. Explicitly, the potential deepens to 

V (x, s = 4) Vi (x) = (ln 4) 2 x (x - 1) (-vr + arcsin y/x) 2 , (39) 

with the particle's speed changing accordingly as a function of x. The return velocity profile 
is now negative, and given by 

v\ (x) = (In 4) y/x (1 — x) (—7i + arcsin y/x) . (40) 
The arcsin in this last expression, as well as in V\, is understood to be the principal value. 
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Moving in this modified negative potential, it now takes the zero-energy particle a finite 
amount of time to travel from x = 1 down to x — 0, as given by At\ = 1. Upon reaching 
x = 0, the exact solution (|27|) exhibits another sudden transition that keeps E = 0, but 
again alters the potential for the return trip towards x = 1. The potential becomes 



VAx) 



V, 



x 



(In 4) ; 



x [x 



1)( 



7T 



arcsm 



(41) 



with corresponding changes in the velocity profile. Again, a finite amount of time is needed 
for the particle to go from x = to x = 1 in the potential V2, as given by At 2 = In (|) /In 2. 

Upon reaching x = 1 again, the switchback process continues. The total energy remains 
at zero, but the potential for the second return trip further deepens, 



V 2 (x) 



V s (x) 



(ln4) s 



x (x 



- 1) (-27T 



arcsm 



(42) 



etc. In order for the particle to follow the interpolating trajectory fl27j) specified by the 
Schroder equation for the s = 4 case, it has to move under the influence of successively 
deepening potentials. 

At later times, the effective potential seen by the particle on its zigzag path between the 
turning points will depend on the total number of previous encounters with those turning 
points. In this sense, the particle remembers its past. Let P be the total number of turning 
points previously encountered on the particle's trajectory. The motion of the particle before 
encountering the next, (P + l)st, turning point is then completely determined, for s = 4, 
by the effective potential (again, arcsin here is understood to be the principal value) 



V P (x) = (In 4) x(x- 1) (-) 



P 



7r + arcsm v x 



(43) 



where [• • • J is again the floor function. That is to say, the potential deepens as P increases. 
The corresponding velocity profile speeds up: 



v P 



(x) = (ln4) (!-*)((-) 



P 



it + arcsm Jx 



(44) 



The particle will either be traveling to the left, with vp (x) < for odd P, or traveling to the 
right, with v p (x) > for even P, with its speed increasing with P. As mentioned previously, 
this effect is clearly seen upon viewing numerical animations of the s = 4 trajectories 0. It 
is also instructive to plot E (t) = v 2 (x (t)) + Vp (x (£)) versus t for various P, with t > tp on , 
to check E = 0, as well as to see how the energy would not be conserved if the potentials 
were not switched. See Figure 1. 

The time for the zero-energy particle to traverse the complete unit interval in x while 
moving through the Vp potential (the Pth "pseudo-half-cycle") is always finite, for P^O: 



At 1 



f 1 dx 


1 








Jo v P (x) 


" hT2 





1 



In 2 



In 



P+l 



P 



This transit time decreases monotonically as P increases, with Atp ~ 

P— »oo 



(45) 

pj^. Starting 



from an initial x, with initial v > 0, the times at which changes in the potential occur, i.e. 
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FIG. 1: E (t) for initial x = 1/2, using potentials V P , with P = 0, 1, 2, 3, & 4. 



the times at which the particle encounters turning points, are obtained by summing these 
transit times and then adding At (x). The transit time sums are simple enough : 



p-i 

N=l 



N 



InP 
Tn~2 



(46) 



Thus, a particle beginning at x, with initial v > 0, will encounter its Pth turning point, 
and the potential Vp will switch on, at time t P on (x) = Ato (a?) + The potential Vp 

will remain in effect for a time span equal to the particle's transit time, Atp = tp + i on (x) — 
tp on (x), i.e. until the next potential Vp+i switches on. 

The potential index P is in fact related to the action for the pseudo-half- cycle [lij • For 
even P, we have / oc ■ ■ ■ , while for odd P, we have / oc J® ■ • • . Thus (in the units 
chosen, the mass is m = 2) 



2 (-If f 
Jo 



vp (x) dx 



1 + P 



+ (-D 



(47) 



where for the first pass through the complete unit interval the action is 

Iq = 2 v (x) dx — (4 In 2) / y/x^/l — x arcsin \/xdx = - (In 2) 
Jo Jo 8 



7T 



(48) 



Moreover, 2 (— v 2 k+i (x) + V2k+2 (x)) dx = 8 (k + 1) Io for integer k, so after the initial 
crossing, the (k + l)st zigzag 1 — > — >- 1 ("pseudo-full-cycle") contributes to the action: 



8(fc+l)I . 



(49) 



The progressive deepening of the potential is also evident in the expanding phase-space 
trajectories for the evolving particle, where the switchbacks enhance the velocity after each 
encounter with a turning point (v — 0). For example, see Figure 2. This Figure provides 
a good opportunity to clarify the differences between conventional Hamiltonian and quasi- 
Hamiltonian systems. The apparent failure to be a conventional Hamiltonian system is due 
entirely to the fact that a given phase-space trajectory seems to repeatedly cross itself at two 
points (x (t) = 0, and x (t) — 1, as shown in the Figure) and at various times. Thus a vector 
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FIG. 2: Phase-space trajectory (x (i) ,v(t)) for s = 4, with x (0) = 1/2 and v (0) = |7rln2. The 
solid portion of the curve is for < t < 3, while the dashed extension is for 3 < t < 3.5. The 
dotted portion of the curve is for negative times, — oo < t < 0, with x (— oo) = v (— oo) = 0. 

field on the phase space needed to describe the motion in the conventional Hamiltonian 
formalism would appear to be singular at those two points. 

These phase-space singularities are tied to the branch point and Riemann surface sheet 
structure of the relevant analytic functions. In fact, an alternate way to visualize the motion 
is as a trajectory on the sheets of a Riemann surface. Consider the particle moving on the 
complex x plane, and not just the real line segment [0, 1], the endpoints of which are now 
branchpoints. There are cuts from +1 to +oo and from to — oo. The particle first 
moves along the real axis, approaches the cut at +1, and then goes around the branchpoint 
such that \fx — > \/x, a/1 — x — > —yl — x, and arcsin — >■ n — arcsin. The particle returns 
along the real axis to the origin and encircles the branchpoint at 0, such that y^x — > —\/x, 
\/l — x — > a/1 — x, and arcsin \fx — > — arcsin yfx. The particle then goes back to +1 and 
goes around it once more such that again yfx —> y/x, a/1 — x — > —a/1 — x , and arcsin — > 
7r — arcsin. The trajectory continues in this way, nipping signs and adding 7rs according to 
the formulae flUD and dH 

For particles with initial v > 0, a third way to picture the dynamics is in terms of the total 
distance traveled by the particle. In this point of view, the successive Vp patch together 
to form a continuous potential V (X) on the real half-line, X > 0, as shown in Figure 3. 
Indeed, the half-line may be thought of as a covering manifold of the unit interval in x, 
with the previous multi- valued functions of x now single- valued functions of X. To avoid 
stagnation at any of the cusps of V (X), either physical or conceptual, in this picture the 
E = right-moving trajectories X (t) may be considered as limits of positive energy right- 
movers, X (t) = lim^o^s (t)- This limiting process corresponds to the encirclement of the 
branch points in the Riemann surface picture. Also, an "extended" phase space trajectory 
involving X (t) and dX (t) / dt gives some improvement over the situation depicted in Figure 
2, with the trajectory no longer self-intersecting. However, the corresponding phase-space 
vector field would still be singular with cusps appearing in dX (t) /dt periodically in X (t), 
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FIG. 3: V(X) where the total distance traveled by the particle is X — x (t = 0) . 

similar to the cusps in the potential V (X). It is again possible to avoid the cusps by 
considering the extended phase space trajectories as limits of positive energy trajectories, 
but we leave that for the interested reader to pursue. 

Finally, in terms of the coordinate transformation introduced in j2[, namely, the angle 



the motion is completely unraveled, at least for s = 4. In that case E = implies 



in a repulsive quadratic potential (inverted SHO). This is the easiest example solvable in 
closed- form by the Schroder method [lj. The exponentially growing solution, 9 (t) = 2 l 9 (0), 
whose values in the real-line cover of the circle are all physically distinct, immediately leads 
to x (t) for the particle moving through the various Vp (x), as is evident from fl27|) . 

To amplify our discussion of the physics encoded in the switchback potentials, we note 
that the chaotic behavior of the s = 4 motion (i.e. the sensitive dependence on initial 
conditions) follows entirely from the behavior of the inverse Schroder function. The sine 
function in (T2"Tj) discards any integer multiple of 2tt. So, for large and increasing times the 
trajectory x (t) | s=4 only depends on digits farther to the right of the decimal in the numerical 
value of the initial arcsin yfx, hence it is sensitive to those initial conditions. Of course, this 
latter fact is well-known in the context of the discrete s = 4 map, and has been discussed 
time and again in the literature for the case where t is an integer (for example, see (^]). 
As is evident from (1271) the discussion carries over mutatis mutandis to the continuous time 
case. The only point we add here is to emphasize that it is which is responsible for the 
effect. Insofar as follows from \1/ and the latter follows directly from the potential (as 
well as the other way round, as discussed in our Introduction) the chaotic behavior follows 
directly from the specific multi- valued form of the potential. 

For other values of s, from ([7]) and we again have x(t) = (s*\l/ (x, s) , s), so 
the distinction between regular and chaotic behavior is again determined entirely by the 



9 = arcsin \fx , 



(50) 




— (In2) 2 # 2 it) = 0, and the trajectory for all t is just that of a particle moving 
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asymptotic behavior of hence the asymptotic branch structure of and hence the 
form of the switchback potentials. The link between the analytic potential function and 
chaotic behavior is arguably somewhat subtle, inasmuch as it requires some computation 
to determine \I/ and but the link is a direct one. Having said that, there remains a 
case-by-case need for inspection of to determine chaos or its absence. At this stage 
and in this paper, we do not purport to produce a universal analysis of all models — we 
only present methods that may facilitate such an analysis. 



4. CONCLUSION 

In summary, in this paper we 

- Use a functional equation to convert nonlinear map equations into dynamical systems 
in the continuous time domain, thereby facilitating the study of such systems. 

- Illustrate this holographic continuum interpolation technique with the well-known 
model of the logistic map, explicitly for the chaotic case s = 4, an example which foreshadows 
an analysis for generic s. 

- Construct an intricate sequence of switchback Hamiltonians which drive the system 
to just those phase-space trajectories resulting from such an analytic interpolation of the 
logistic map, and produce the familiar chaotic behavior for s = 4. 

For nontrivial systems such as the chaotic s = 4 logistic map illustrated here, in order to 
produce the global features of the trajectories the potential of the underlying Hamiltonian 
system must change at each switchback of the motion, as occurs when turning points are 
reached in finite time. Therefore we have called these systems quasi- Hamiltonian. The 
trajectories dictate and completely determine the succession of switchback potentials, in a 
richly extended analogy of inverse scattering techniques. 

For values of s between 2 and 4, a numerical analysis pursuant to the techniques in 
this paper leads to quasi-Hamiltonian systems involving switchbacks among sequences of 
potentials, eventually producing trajectories that move between elaborate sets of turning 



points, in general. This analysis is reported in full in a companion paper 16]. However, 
to make the present paper more self-contained, and perhaps more convincing as to the 
generality of the switchback effect, we have included some details for one other case (s = 3) 
in an Appendix. 

To conclude, we encourage experimenters to look in various settings for the particular 
type of continuous interpolates that we have described, including fluid dynamical experi- 
ments, ecological systems, and, more generally, any situation where an inherently continuous 
parameter (time, scale, etc.) has been routinely sampled only at predetermined discrete val- 
ues. 
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FIG. 4: R\£, (n) — s/4 for s = 3, plotted for 2 < n < 400, to show convergence towards = 3/4. 



APPENDIX: SOME NUMERICS FOR THE s = 3 CASE. 



The coefficients in the series expansion of \1/ may be computed recursively for any s, from 
( I2T1) . While not immediately recognizeable as the series expansion of any well-known func- 
tion, nonetheless, the coefficients may be used to estimate the radius of convergence of the 
series (|20|) by the simple ratio test: Rq, = lim^oo (n), with R^ (n) = \d (n — 1) jd (n)\. 
To illustrate this consider the case s = 3. For this case numerical application of the ratio test 
leads to the data in Figure 4. From the numerical data we infer that Rq = 3/4 = s/4| s=3 , as 
indicated earlier in (1221) . While this numerical demonstration is not a proof of convergence, 
of course, we believe it provides compelling evidence for convergence. 

Other values of s produce numerical data qualitatively the same as the s = 3 case, with 
only one noteworthy difference: While for s > 1, the convergence of the numerical data to 
the expressions in f[2"2l is from above, for s < 1 the convergence is from below. 

A similar numerical calculation for the coefficients in the expansion of the inverse function 
ffTB"]) suggests (x) is an entire function of x, for s = 3, much like the closed-form results 
in f|T7|) . That is to say, R^-i = 00 for s = 3. 

Combining the series data for ty' 1 , and V with the functional equations obeyed by 
them leads to an accurate numerical determination of the functions and the switchback 
potentials, for any s. The methods are described in detail in 161 ] . Numerical results for 

are shown in Figure 5 along with the expected asymptote, 
^l/" 1 = 2/3. This interval in t corresponds to plotting (x) for < x < 10 4 . 

Numerical results for the first few potentials in the sequence, for s = 3, are shown in 
Figure 6. For this particular case, the sequence of potentials does not deepen progressively, 
but rather the sequence converges, albeit very slowly, onto the stable fixed point of the 
map, namely, = 1 — ^\ s _ 3 = 2/3, shown in the Figure as a small circle on the x-axis. 
Hence the zero-energy particle trajectory will also converge onto this fixed point, rigorously 
conserving E = in the process. 

Although the particle trajectory is not chaotic for s = 3, nevertheless the potentials are 
somewhat more exotic than the s = 4 case in the sense that there are an infinite number 
of branch points for the s = 3 case. These may be obtained by iterating the action of 
the map 0111 = 1/2. That is to say, \ -> | = Jr ^ = J -» iff = 3 3 ^ §f§ = 
— > ■ ■ ■ . These branch points are shown as the turning points of the potential 



047 67 
' 2 16 



sequence, 0.75, 0.562 5, 0.738 281, 0.579 666 



represented by colored dots on the x- 
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FIG. 5: tf -1 (i s*) | s=3 , in blue, for -4 < i < 9, from series plus functional methods. 

axis in Figure 6. The branch points may also be seen by flipping the graph of in Figure 
5 about the SW-NE diagonal, in the usual way, to exhibit the various branches of 

This complicated branch structure for s = 3 extends to other cases with 2 < s < 4 
and most likely accounts not only for the failure to render them in terms of simple, known 
functions, but also for some of the associated wide-spread mystique that surrounds \l/ and 
^f -1 in these cases. Of course, s = 2 and s = 4 only have branch points at x = 1/2, and at 
x = and 1, respectively, as evident in the closed-form results ( TT7|) and as discussed in the 
text, so these cases are relatively bland by comparison. 
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